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INTRODUCTION 



Clusters of galaxies have long been recognized as remark- 
able laboratories to test cosmological theories. They are 
the largest gravitationally bound objects in the universe, 
thought to have formed from the eventual collapse of ini- 
tially overdense dark matter perturbations. Their abun- 
dance and large scale properties are sensitive to the ex- 
pansion and initial conditions of the universe, making 
them excellent tools to constrain cosmological models. 
On smaller scales, the physics of clusters is dominated by 
complex baryonic processes such as gas cooling, star for- 
mation, and feedback from supernovae and active galactic 
nuclei. In order to use clusters as standard probes of the 
geometry and dynamics of the universe, it is necessary to 
reliably model these processes and distinguish amongst 



ABSTRACT 

We present a statistically-optimal and model-independent method to extract the 
pressure profile of hot gas in the intracluster medium (ICM). Using the thermal 
Sunyaev-Zeldovich effect, we constrain the mean pressure profile of the ICM 
by appropriately considering all primary cosmic microwave background (CMB) 
and instrumental noise correlations, while using the maximum resolution and 
sensitivity of all frequency channels. As a first application, we analyze CMB maps 
of WMAP 9-year data through a study of the Met a- Catalogue of X-ray detected 
Clusters of galaxies (MCXC). We constrain the universal pressure profile out 
to 4i?5oo with I5(T confidence, though our measurements are only significant 
out to -R200- Using a temperature profile constrained from X-ray observations, 
we measure the mean gas mass fraction out to i?2oo- Within statistical and 
systematic uncertainties, our constraints are compatible with the cosmic baryon 
fraction and the expected gas fraction in halos. While Planck multi-frequency 
CMB data are expected to reduce statistical uncertainties by a factor of ^ 20, 
we argue that systematic errors in determining mass of clusters dominate the 
uncertainty in gas mass fraction measurements at the level of ^ 20 percent. 

the ICM. rj Historically, most of our observational un- 
derstanding of the ICM has come from X-ray observa- 
tions, leading to a fairly consistent picture of the scaling 
and structural properties of low-redshift clusters (see e.g. 
^Bohringcr ct al.'2007: V ikhhnin et al.'2006: Arnaud ct al. 
|2010,,Pratt et al.,2007,,Sun'et al. [2011,, Pratt et al.,2009f7 
The thermal Sunyaev-Zeldovich (tSZ) effect ( Sun- 
|yaev fc Zeldovich||1972[ ) is another important probe of 
the ICM: as the Cosmic Microwave Background (CMB) 
photons inverse-Compton scatter off of the hot electrons 
in the ICM, their blackbody spectrum is distorted. The 
tSZ effect has the unique property that its signal is inde- 
pendent of redshift, making it a powerful observational 
tool for detecting clusters at cosmological redshifts, and 
hence a promising cosmological probe of dark energy 
(e-g 



Birkinshaw 1999 Carlstrom et al. 20021. Within 



different feedback mechanisms (see Kravtsov & Borgani 
2012[ [Allen et al.|20li| for recent reviews). 



the last few years, cluster surveys exploiting the SZ effect 
have started delivering cluster samples (e.g., Staniszewski 
et al.|2009| [Marriage et al.l2011[|Wilhamson et al.|2011 



The main baryonic budget of clusters is a hot 
plasma of ionized hydrogen and helium in the intracluster 
medium (ICM), making it the natural target for study- 
ing the complex astrophysical processes at play. This 
virialized plasma emits bremsstrahlung radiation in X- 
ray, making it possible to probe the dense regions of 



Planck Collaboration VIII[[2011a l as well as constraints 
on cosmological parameters (Vanderlinde et al.|2010 Se- 



hgal et al. 20111. Another important feature of the tSZ 



^ This is because X-ray luminosity is proportional to the 
square of gas density (e.g. [Weinberg|2008[l . 
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efTect is that it is directly proportional to the integrated 
pressure of free electrons along the line of sight, which 
makes it a powerful probe of the ICM in the outskirts 
{r > R200), where X-ray emission is extremely faint. 

Resolving the tSZ signal for individual clusters re- 
quires high resolution CMB measurements, which have 
become available only in the recent years. In fact, the 
only all-sky CMB survey with high enough resolution 
and sensitivity to detect individual SZ clusters is Planck. 
Even with Planck's sensitivity, it is necessary to combine 
the tSZ signal from many clusters to meaningfully con- 
strain physical quantities of interest, such as baryonic 
mass fraction (Planck Collaboration V 20131. Luckily, 



there is fairly strong evidence from X-ray observations 
and numerical simulations that clusters are self-similar 



to a good approximation (see e.g. Nagai et al.|2007 Ar- 
|naud et al.||2010| . This fact justifies combining SZ sig- 
natures of many clusters to obtain constraints on the 
mean ICM properties. There have been quite a few ef- 
forts in this direction over the past few years. By analyz- 
ing WMAP-1 (WMAP-3 respectively) CMB data for 116 



(193 respectively) X-ray detected clusters, Afshordi et al. 
( |2005[ |2007[ ) (respectively) provided constraints on the 
ICM pressure profile out to ~ -R20o- Other similar works 
include WMAP-3 stacking of over 700 clusters by [Atrio^ 
[Barandela eTaTlpOOSl), WM AP-7 analysis of 175 Planck 
ESZ clusters ( Ma et al. 2013 1, and SZ measurements of 15 



massive X-ray selected clusters obtained with the South 
Pole Telescope ([Plagge et al.|2010[ ). Most notably, [Planck| 
Collaboration V ( 2013 1 have studied the tSZ signal of 62 



low-redshift massive clusters by using CMB data from the 
Planck satelite, constraining the mean pressure profile of 
the ICM out to 3-R500 with unprecedented precision. 

The practice of averaging signals from many clus- 
ters goes under the title "stacking" . The basic idea is the 
following: the main sources of uncertainty in extracting 
the tSZ signal are the primary CMB anisotropies and 
instrumental noise. Since these sources of noise are ran- 
dom in nature, they "drop out" if the temperature profile 
around many clusters is averaged over. This procedure, 
in its various forms, is not statistically optimal for mul- 
tiple reasons. Firstly, large-angle correlations of primary 
CMB fluctuations are ignored when stacking. Secondly, 
it is not clear how contributions from different clusters 
should be optimally weighed in the averaging process. 
Typically, different weighing methods are adopted to see 
whether the effect on the extracted tSZ signal is signif- 
icant or not (e.g. Atrio-Barandela et al.| [2008). Thirdly, 
when using multiple frequency channels, the final resolu- 
tion of the reconstructed tSZ map is determined by the 
lowest resolution of the combined frequency maps (e.g. 



Planck Collaboration V 20131. Finally, the 3D pressure 
profiles are usually obtained a posteriori by deprojecting 
the tSZ signal, which may lead to noise amplification. 
For these reasons, stacking procedures either result in an 
underestimation of error or loss of statistical information. 
We believe all the aforementioned shortcomings of 
stacking procedures can be overcome with the method- 
ology we have formulated in this paper, which is more 
in line with optimized template fitting procedure of |Ko-| 
matsu et al. (20111. Our analysis includes an all-sky 



which appropriately takes into account primary CMB 
and noise correlations on all scales, while using the maxi- 
mum resolution and sensitivity of all channels to their full 



potential. Furthermore, following Afshordi et al. (20071 



our method is completely model- independent, thus elim- 
inating any systematic uncertainty associated with the- 
oretical modelling of the ICM. In this paper, we apply 
our formalism to WMAP-9 CMB data using the Meta- 



Catalogue of X-ray detected Clusters of galaxies (Pif- 
faretti et al.|201l K In a companion paper, we will repeat 
our measurements using CMB data from the Planck mis- 



Our paper is organized as follows. In Sections [271] and 
|2.2| we review the tSZ effect and the concept of a univer- 
sal pressure profile, describing how model-independence 
can be achieved. Sections 12.31 and 12.41 contain our main 
statistical and numerical methodology, outlining in detail 
how the mean pressure profile of ICM can be optimally 
constrained. Section |3] describes the CMB data and clus- 
ter sample we use to test our methodology, and is followed 
by a discussion of the resulting pressure profiles in Section 
|4.1[ Section [4. 2| presents gas mass fraction measurements 
of various subsamples of our cluster catalogue. We dis- 
cuss future work and how we anticipate our results to 
improve by using Planck CMB data in Section [5] before 
concluding our findings in Section [6] 

Throughout this paper, we assume a ACDM cosmol- 
ogy with present matter density fim = 0.3, dark en- 
ergy density Ha = 0.7, and Hubble parameter Ho = 
100 h km s'^Mpc'Vith h = j^hro = 0.7. We also de- 
note the normalized Hubble parameter at redshift z by 



+ nA. 



2 EXTRACTING THE PRESSURE PROFILE 

This section contains the statistical and numerical 
methodology we use to extract the mean ICM pressure 
profile from a full-sky CMB experiment. Section [2T] re- 
views the tSZ effect and how it is related to the electron 
pressure profile. In Section [2.2[ we reduce the problem of 
finding the exact profile of each cluster to a single, univer- 
sal up to normalization, pressure profile. In Section [2. 3| 
we derive the maximum likelihood estimator of the pro- 
file and its covariance matrix. Finally, in Section [2. 4| we 
describe how the components of the estimator are in prac- 
tice computed numerically. 



2.1 tSZ Effect Model 

The contribution of the tSZ effect to the CMB temper- 
ature anisotropy at frequency i^ and location n on the 
sky is proportional to the integral of the electron pres- 



sure along the line of sight: (see §2.5 of Weinberg (20081 
for a derivation) 



hu 

ksTcMB 

F{x) = a;coth(2:/2) -4, 



rUcC^ 



I 



aifiic \''fi) ^ 



multi-channel fit to the mean pressure profile of the ICM 



where ctt is the Thomson scattering cross-section, rrig 
is the mass of the electron, c is the speed of light, ks 
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is the Boltzmann constant, h is the Planck con stant, p 
TcMB = 2.725 K is the mean CMB temperature ( [Mather 
et al.|[l999l, and Pc{ln) is the pressure of free electrons 



along the line of sight direction n. 

Our task is to constrain P^ through the tSZ effect. 
We will assume that Pe is spherically symmetric to a good 
approximation and denote the pressure profile of the a**^ 
cluster by P^ (r). Furthermore, since it is not possible 
to constrain a continuous function without introducing 
model-dependence, we consider spherical bins around the 
centre of clusters, in each of which the pressure is assumed 
to be constant: 



(Pi 



(a) 



pr'ir) = { 



3(a) 



,{") 



if < r < r^°' 
if 




(a) 



<r< 4"^ 



(2) 



if r 



(a) 



< r < r 



(a) 



Figure 1. The fc*'^ bin of the a*'' cluster. The contribution 
of this bin to the temperature anisotropy of the CMB at fre- 
quency u and location n on the sky is given by P^."^ tj^ ("■: i^), 



Here P^"^ , . . . , Pjl^' are all constants with units of pres- ^^ere 4"^ (fi; u) is defined in Equation ^ 



sure, r is the radius away from the centre of cluster, and 
A'';, is the total number of bins. The value of the pressure 
in each bin may be better understood as the volume- 
weighed average of the pressure in that bin. With these 
simplifications. Equation ([l]) may be written as 



fc=l a=l 



(3) 



where N^ is the total number of clusters and t 
given by: 



(")^v^. 



4">(n;^) 



o^TTc^ 



•-F 



hv 



X 2 



;{o) 



ksTcMS , 
if K;('")(n) <r\^^ 
if ri"' <»('') (")<^Si (4) 

The functions used in Equation (Mk are defined as follows: 



;ri(n)-/r(n) 



l^\n)^^[ri^\^-[w^-K'h)\\ 


(5a) 


w''')(n) = d("'sin(6l<"'(n)), 


(5b) 


cos(6)'">(n)) = n-ri<"'. 


(5c) 



where d'"' is the angular diameter distance to the a"^ 
cluster and n'°' is the unit vector pointing to its centrerj 
Fig. [l]shows the basic geometry that underlies Equa- 
tions |3][5|. 

2.2 Universal Pressure Profile 

In principle, the analysis that will follow can be used to 



optimally estimate all parameters P{ 



(a) 



,P^;>. Unfor- 



tunately, this is too computationally-intensive for a large 



sample of clusters, given the large angle correlations of 
primary CMB anisotropics and the current resolution of 
CMB experiments. However, there is fairly concrete evi- 
dence that the pressure profile of the hot gas in clusters 



is self-similar (see e.g. Nagai et al. 2007 Arnaud et al. 
2010 1. This means that for a given cluster, there is a self- 
similarity scale re such that the pressure profile takes 
the form 



Pr(r)=P(r/rr)^c^ 



(a)\ p(a) 



(6) 



where P^ is a constant characteristic pressure of the 
a* cluster, and P(a;) is the so-called universal pressure 
profile. Within this context, it is natural to construct the 



M) 



radial bins so that r^"'' = Ukr)^' , where {nk}k=\~Ni, are 



M) 



positive numbers satisfying n\ < n2 < ■ ■ ■ < tim^ ■ With 
these assumptions, discretization of the electron pressure 
profile amounts to Pj^ = PfcPc , and the tSZ contribu- 
tion (|3| to the CMB anisotropy takes the form 



&T''^{h-u) = Y.¥kt'^^\n), 



fe=i 



f('")r^ 



(n) = ^pW4'">(n;^). 



(7a) 



(7b) 



If the characteristic scales r^ and Pc are fixed by 
external observations (such as X-ray's) for each clus- 
ter, our task is simplified to finding best estimate values 
(and their associated uncertainties) for A*';, parameters: 
Pi,...,Piv,. 

We use P5Q0 as the self-similarity length scale of the 



M) 



?(<>)• 



(a) 



a cluster (i.e. we set re = Psoo)- The quantity R)^ 
defined as the radius up to which the matter density is 
A times the critical mass-density of the universe: 



^ We've used h to denote Planck's constant only in Equations 
^\ and (J4|. Throughout the rest of our paper, h is the reduced 
Hubble constant. 

^ Angular diameter distance to the a*'^ cluster with redshift 
z is given by d(")(z) = ^/^ L^ ''^ 
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M 



(a) 



Pm'(r-)47rr dr 



= A X -ARlV X Pcr^t{Za) 



(8) 



where Za is the redshift of the a*'^ cluster, and pcrit{z) 
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„M 



We consider 8 bins with radii rj, — 0.5k x -R5001 



?(<>) 



3-H"(z) 
SttG 

where k £ {1,2,..., 8}. This is equivalent to setting nt 
0.5k. We use two proposals for the characteristic cluster- 
dependent pressure Pc ■ 



,(") 



1.65 X 10~^E{za) 



8/3 



hn 



M. 



(a) 



3 X lOi^/i-giMg 



2/3+s 



keV cm 



(9) 



with 5 — and S = 0.12. The former corresponds to 
the mass variation expected in the standard self-similar 
model based purely on gravitation (see Nagai et al.|2007 



Arnaud et al.|[2010|. The S — 0.12 case is a modification 



to the standard self-similar model proposed by lArnaudI 



et al. ( 2010 1, which is a first approximation to quantifying 
how the mass scaling of Pc changes with radius in the 



REXCESS sample ( [Bohringer et al.|2007[ ). Using the 6 ■■ 
0.12 scaling makes for a meaningful comparison of our 
results with those of [Planck Collaboration V| ( |2013[ ), since 
this is what is used in their analysis. 



2.3 Statistical Methods 

We use the principle of maximum likelihood to find best- 
estimate values for the parameters Pi , . . . , Pjv^ . (For tem- 
plate fitting on CMB sky, see e.g.|Gorski et al.|1996l|Jaffe| 



et al. 2004 Komatsu et al. 20111. We assume that the 



only contributions to the temperature anisotropies of the 
CMB are the primordial anisotropies 5T^ , the tSZ effect 
5T (given by Equation (|7|), possible residual monopole 
and dipole components, and the instrumental noise. Fur- 
thermore, we assume that primordial anisotropies follow 
the statistics of an isotropic Gaussian random field, for 
which we know the angular power spectrum C^.n 

Let Lp {Lu) denote the set of all pixels (frequency 



channels) we wish to use in our analysis. In Appendix Al 
we show that the log-likelihood of the observed temper- 
ature fluctuations {^T'ii^li'g^", given the binned pressure 
profile {Pfe}, is analytic and equal to 



\x\{5T.A'iiL;\{^k}) ^ 



hsT-ST^^fC-^ST-ST^'^), (10) 



where 



C = Cs + Cf 



(11) 



Here Cs is the covariance matrix of the primary CMB 
fluctuations and Cn is the covariance of the instrumental 
noise. The former is related to the angular power spec- 
trum d through 

[Cs]i^,iV = {STi^ STi,^,) = 

E ( -^ ) aBi,B,,,w!Pe{fi, ■ n,,), (12) 



* We use the CAMB code to generate the expected values 
of Ci's for the WMAP concordance ACDM cosmology | |Lewis| 
let al. 120001 IHowIett et al 



|2012[|. 



where Biv is the spherically averaged beam transfer func- 
tion for the mode £ and frequency channel i/, We is the 
spherically averaged pixel transfer function, fii is the sky 
direction of the i* pixel, and Pe is the ^-th degree Leg- 
endre polynomial. We use Imax ~ 2 x N^ide throughout 
our analysis, where A^sWe = 512 is set by the HEALPix 
( Gorski et "aL]|2005[ ) resolution of the WMAP sky maps. 
In the case of WMAP, the instrumental noise is 
largely uncorrelated both between pixels and between dif- 
ferent frequency channels. Its properties are adequately 
modelled by a Gaussian distribution with covariance ma- 
trix 

In Equation ( 10 1, we have used STfJ' to stand for a pix- 
elized version of Equation ([7| : 



5Tl^^Y.^utt^{i), 



(14a) 



4^^W = E E {t^:\raBuMlYlr^{n^), (14b) 



i = m = -i 



f(■')^ 



where (ij. )im are the spherical harmonic coefBcients 
of t^ {fi). In practice, we generate the templates at 
HEALPix resolution 12 (iVsido = 4096), and then down- 
grade to resolution 9 (A'^sidc ~ 512) P] We then convolve all 
templates with instrumental beams to obtain the quan- 
tities t''^\i). 

Let us note again that P^ is the universal pressure 
of the fe**^ bin and t^^ {i) is the coefficient that multi- 
plies it at pixel i and frequency channel v. Assuming a 
uniform prior on (Pi, . . . ,PjVi,), the posterior probability 
function of these variables is also a Gaussian distribu- 
tion which peaks at the maximum of the log-likelihood 



function given in Equation (10 1, which is 

Pfc = Ei^^'u.fe'/^fe'' 



Q!fc,fc' 



Pk 



E E4^'»4'"(^ 



(15a) 
(15b) 
(15c) 



u^L^ i£Lp 



where we have introduced the inverse covariance weighed 
template 



xi^^ii)- E E [c"']-.'^'4"^''(0 



(16) 



u' ^Lfy i' ^Lp 



Finally, the covariance matrix Cp of {Pfe} is determined 
by the Hessian of the log-likelihood ( |10| |. Its matrix ele- 
ments are 

--[cx-\^,,. (17) 



^ To be more specific, following the notation of Section |2.2[ we 
give the value V^ (fii) to the i-th pixel of the fc-th template 
in frequency channel v, where fn is the centre of the i-th 
pixel at a resolution 12. (At this resolution, the first radial bin 
of all MCXC clusters occupies at least one pixel.) We then 
downgrade the templates to resolution 9 and finally convolve 
all templates with the instrumental beams. 
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This matrix does not involve further computations as it 
is already required to obtain the best estimates {Pfc}. All 
the measurements on the pressure profiles and their at- 
tached error bars are obtained using only the expressions 
indicated in this Section. 

In Appendix |A4[ we have shown how any residual 
monopole or dipole contribution can be conveniently ac- 
counted for in this formalism. 



2.4 Numerical Methods 



It is clear from Equation ( 15 1 that all quantities of inter- 



est can be calculated once the weighed templates Xj;,'^ 
are known. This is impossible to achieve by direct com- 
putation, which would involve inversion of the full covari- 
ance matrix C. In this section, we will describe how we 
compute Xj^'^ numerically. 

We start by establishing some notation. The set of 
all masked (unmasked) pixels is denoted by Lp (Lp), so 
that L = LpULp contains all pixels in the sky. We denote 
the total number of pixels (i.e. the size of L) by Nt- 

In the statistical modelling described in Section [2. 3| 
the covariance matrices Cs and Cn are only computed 
on the observed pixels Lp. However, numerical manip- 
ulation of these matrices is more efficient in harmonic 
space, using the spherical harmonic transform, which it- 
self requires the knowledge of all pixels. Therefore, it is 
advantageous to compute the quantities of interest by ex- 
tending the domain of Cs and Cn to the entire sky. We 
refer the reader to Appendix |A2| for details on how this 
can be achieved and state the final result here (see also 
[Wandelt et al.|2004|[Lavaux et al^[20T2l ). 

We encode information about masking of pixels into 
a diagonal Nt x Nt matrix M, with elements Mu = if 
i G Lp, and Mu — 1 ii i £ Lp. Furthermore, we let S be 
the full pixel to pixel covariance matrix due to primary 
CMB fluctuations: (i, j G L) 



[S]., = I] C7, f '^^-^\ Piin, ■ h. 



47r 



(18) 



N^ the pixel-to-pixel covariance matrix of the instrumen- 
tal noise: 



[.LN[/Jjj — IT'ii^Oi^ 



(19) 



and B^ a model of the complete beam (pixelization and 
instrumental): 



[B^]ij — ylpix y ^ 



21+1 
4tt 



Be,WiPe{n,-nj), (20) 



where Apix is the area of one pixel (which is equal to 
A-k/Nt for all pixels in the HEALPix scheme). We further 
define the square root of S: 



[S ' \ij = \/Apix 2_^ ^cv 



21+1 
4tv 



Pi,{h,-nj). (21) 



variance matrices, temperature data, and templates: 

M 



¥u=Y,{a.-^)k,k'Py, (22a) 

"*.'='= EE*i'''»-^^"''»' (22b) 

h=Y.Y.^T,''^i^\i). (22c) 



i^&Lu ieL 



where 



4""w= E E[G^-']-'*i""(0. 



(23a) 



v'eL^ i'eL 



mn;;^mb^s^''^d"^s^''^b^'Mn;/m, (23b) 

D = 1 + S^''^ [ Y^ Bj.MN-^MB^ ] S^/^ (23c) 
The computation of Xl!^ is now reduced to solv- 



es) 



D-if^''\ where f^"^ 



ing for the quantities gj, 
S^/^B^MN^^Mt^''^ The other operations may be done 
trivially as all involved operators are either diagonal in 
pixel space or in harmonic space. Computing gj:" is 
equivalent to solving the equation Dg^" = f J." , for which 
a number of numerical techniques are available. We use 
the algorithm of the conjugate gradient method with pre- 
conditioning (see e.g. Shewchrik||1994 1, which is an iter- 
ative prescription for solving large linear systems equa- 
tions of the type 



Ace = y. 



(24) 



To speed up the convergence, we construct a pre- 
conditioner matrix Hq^ (essentially an approximation of 
D~^) as follows: the block corresponding to all harmonic 
modes with ^ < 60 is taken to be the exact inversion of 
the same block in D, which is computed using a Cholesky 
decomposition. The rest of D^^ is taken to be diagonal, 
the elements of which are reciprecals of the corresponding 
diagonal element of D. This preconditioner has already 



been used in other works (see e.g. Eriksen et al. 20041 



We stop the conjugate gradient algorithm whenever the 
relative error 



\AXn - y\\2 

^l\y\\2 



(25) 



Equation ( 15 1 may now be rewritten using the full co- 



is less than a specified threshold. In the case of this work, 
we take e„ < 10~^. We have checked that changing this 
threshold to 10~^ does not change the results, indicating 
that the solution has indeed converged (see Appendix [b] 
for detailed convergence tests). 

Finally, we note that all {g^ } are fully indepen- 
dent and thus may be computed in parallel. We fully 
employ this property. Our software, ABYSS (the spher- 
ic Al BaYesian Statistical Sampler), runs in 25 hours and 
52 minutes on seven nodes (16 cores) to solve for the 
12 templates on an Intel Xeon E562Q. We note that the 
monopole and dipole take significant more time to reach 
the same level of precision as the other maps. 
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3 DATA 

3.1 CMB Data 

We use co-added inverse-noise weighted data from nine- 
year maps observed by WMAP at 41 GHZ (Q-band), 
62 GHZ (V-band), and 94 GHz (W-band) p] These maps 
are foreground cleaned ( Bennett et al.|[2012 | and are at 
HEALPix resolution 9 {Nstde = 512). The standard de- 
viation of the pixel noise in each map is given by (using 
notation of Section 2.4 1 



J-) 



vw 



(26) 



where v £ L^ = {Q, V, W}, i G L, and 4*^' = 2.188 mK, 
t7^^' = 3.131 mK, cr^*^' = 6.544 mK. The number of ob- 
servations A'^"'* at pixel i is included in the maps avail- 
able from the LAMBDA website. In all of our analysis, 
we use the 'extended temperature data analysis mask' 
to exclude foreground-contaminated regions of the sky 
from the analysis. The beam transfer function for every 
differencing assembly is also provided on the LAMBDA 
website. For a single value of I, we average beam transfer 
function values for all differencing assemblies belonging 
to the same frequency channel. This is how we obtain the 
quantities Bi^ introduced in Section [2. 3[ 



3.2 Cluster Sample 

We use the Met a- Catalogue of X-ray detected Clusters of 
galaxies (M CXC) to extract the u niversal pressure profile 
of the ICM jPiffaretti et al.|201l[ ).nThe MCXC provides 
(amongst other quantities) sky coordinates, redshift, and 
M500 data for all clusters. With a few exceptions, lumi- 
nosity is used as a mass proxy for all clusters (see equation 
(2) of |Piffaretti et al.|[20Tl| . 

We perform our analysis on all 1743 MCXC clus- 
ters, as well as a subsample of 162 clusters whose first 
radial bin {= 0.5 x i?5oo, as discussed in the next Sec- 
tion) is resolved by the W frequency channel of WMAP. 
More specifically, we obtain this subsample by requir- 
ing d(z)6i'^) < 0.5i?500, where O'^^^ = 0.12° is the ef- 
fective angular radius of the (averaged) W-channel de- 
tector beam. We will refer to clusters in this subsample 
as resolved MCXC clusters. Fig. |2] shows the redshift- 
i?50o distribution of all MCXC clusters, differentiating 
between the resolved and unresolved ones. The redshift 
of aU (resolved) MCXC clusters range from 0.0031 - 1.26 
(0.0031 - 0.077) with a median of 0.14 (0.028), and their 
masses range from jptf^ = 0.0096-22.1 (0.0096-7.27) 
with a median of 1.77 (0.86). 

Since most MCXC clusters cannot be resolved, one 
expects numerical uncertainties to become important. 
This is why we have chosen to study a subsample in which 
all clusters are resolved. However, even unresolved clus- 
ters contribute to the tSZ signal, especially in the outer 
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^ http://lajnbda.gsfc.nasa.gov 

^ All information about MCXC clusters may be found here: 

http : //vizier ■ cf a . harvard ■ edu/vlz-bin/VizieR?- sour ce= | 

'j/A+A/534/A109l 



R,oo (Mpc) 

Figure 2. The redshift-i?5oo distribution of MCXC clusters. 
The red points represent clusters whose first radial bin is re- 
solved by the W frequency channel of WMAP. The three lines 
plot (^(2)6'"' for different WMAP frequency bands, where 

S'"^' = i/r2('^)/7r is the angular radius of the disk with the 
same effective area as the detector beam in frequency channel 
u. Here O'"'' is the beam solid angle of frequency channel u, 
which is provided on the LAMBDA website: OW) = 0.51^, 
a(V) = 0.352, ^^^ q(W) ^ Q 222 (deg^). The curve labeled 
'Pix' is constructed similarly and reflects the resolution asso- 
ciated with pixelization. 



bins. Therefore, the price one pays for ignoring unresolved 
clusters is statistical information. We have analyzed both 
samples to see how this trade-off manifests itself in prac- 
tice. 

We also analyze subsamples of MCXC clusters 
binned according to mass. This allows us to study the de- 
pendence of various quantities, such as pressure and gas 
mass fraction, on the mass of clusters. Table [T] shows the 
mass range and number of clusters in every bin. We have 
subdivided the resolved MCXC clusters into three mass 
bins, and the entire MCXC sample into four bins. These 
bins have been chosen so that they lead to roughly sim- 
ilar signal-to-noise properties, characterized by the null 
chi-squared of pressure measurements. We have excluded 
the 45 most massive clusters because none of them are 
resolved, resulting in a measurement with extremely low 
significance and a nearly degenerate covariance matrix. 



4 RESULTS 

In this section, we describe our two main results. Sec- 
tion 14.11 discusses WMAP constraints on the universal 
pressure profile P, and Section [4.2| includes our gas mass 
fraction measurements. 



4.1 WMAP Constraints on the Universal 
Pressure Profile of the ICM 

Fig. [3] shows the result of our pressure measurements as 
applied to all MCXC clusters, as well as the resolved 
subsample defined in Section |3.2[ As was mentioned in 
Section [2. 2| we perform our analysis using the standard 
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Table 1. Binning MCXC clusters according to their mass. 



Bin number M500 range (IO^^Mq) Number of clusters 



1 


0.0096-2.41 


1140 


2 


2.41-4.17 


364 


3 


4.18-5.31 


124 


4 


5.32-7.27 


70 



(a) All but the 45 most massive MCXC clusters. 


Bin number 


M500 range (IO^'^'Mq) Number of clusters 


1 
2 
3 


0.0096-2.71 138 
2.83-4.56 15 
5.17-7.27 9 



(b) Resolved MCXC clusters. 



-,M 



2/3 



self-similar model Pc oc Mjqq {5 — in Equation ||9| 
as well as the modified scaling Pc oc M^qq ' (<5 



3a 



0.12 in Equation ([9l). These results are shown in Fig 
and Fig. |3b| respectively. There is essentially no signal 
beyond 1.5i?5oo. The best fit pressure values even become 
negative for some bins in this regime. We have decided 
not to impose positivity of pressure as a prior in order 
to keep the statistics Gaussian and not spoil the analytic 



results ( 15 I and ( 17 1. Repeating these measurements with 



Planck CMB data is expected to provide a significantly 
tighter constraint on the universal pressure profile (see 
Section [s] below) . 

In Fig. |3a| we have compared our pressure measure- 
ments with two sets of simulations (Nagai et al._ 2007| 
Battaglia et al.|[2012 1 , which include treatment of radia- 
tive cooling, star formation and energy feedback from su- 
pernova explosions. Battaglia et al. (20121 also account 



for feedback from active galactic nuclei, while [Nagai et alT] 
([2b07|) consider the effect that electrons and ions are not 



kept in thermal equilibrium in the outskirts (Rudd & 
Nagai |p009 1 . Comparison of our measurements with the 



simulated profiles of Nagai et al. ( 2007 1 is straightforward 
because they compute the exact same quantity. Battaglia! 



et al. ( 2012 1, however, use i?200 as the self-similarity scale 
and also consider the variation of P with mass and red- 
shift. In this case, we use the C200 — M200 relation of 



Bhattacharya et al. (2011 1 to compute R200, and use the 



fitting formula of Battaglia et al. (20121 (equations (11- 
12) and table 1) to compute P for all MCXC clusters. 
In Fig. [3a| we have plotted in green the average of these 
profiles, as well as the standard deviation about their 
mean. Where there is signal, our pressure measurements 
are slightly more consistent with those of [Battaglia et aL] 
(20121. Due to the large statistical uncertainties, how- 



ever, we cannot meaningfully discriminate between the 
two. In an upcoming companion paper, we will return to 
this question in greater detail once we apply our method- 
ology to Planck CMB data. 

In Fig. |3b[ we have compared our pressure mea- 



surements with those of Planck Collaboration V (20131 



Given that we use different radial bins, and more impor- 
tantly that Planck's measurements are a lot more pre- 
cise, it suffices to compare our measurements with their 
best fit generalized Navarro-Frenk White (GNFW) pro- 




(a) Standard self-similar scaling of pressure with mass (5 = 
in Equation Kj). The shaded areas mark the dispersion 
about the mean profiles of simulated clusters from |Nagai | 
et al.|||2007|l (gray), and|Battaglia et al.|||2012^ (green). 



hJh hfn WMAP9-AII MCXC Clusters 
hJh hfn WMAPg-Resolved MCXC Clusters 
— Planck Best Fit (PBF) 
• • Binned PBF 




-0.2 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

^'=r/ 1^500 

(b) Modified self-similar scaling of pressure with mass 

(<5 = 0.12 in Equation l|9l). The black curve is the best fit 

GNFW profile to pressure measurements of Planck. The 

black points are obtained by a volume weighed average of 

Planck's best fit profile over our radial bins. 

Figure 3. WMAP-9 constraints on the universal pressure pro- 
file P of the ICM. The blue (red) data points are the resulting 
pressure profiles for all (only resolved) MCXC clusters. A clus- 
ter is considered resolved if its first radial bin subtends a solid 
angle larger than the effective beam area of the W frequency 
channel (see Section [3. 2[ l. Defined in Equation (J9l, S charac- 
terizes deviation from the standard self-similar model. 5 = 
corresponds to the mass variation expected in the standard 
self-similar model (see [Nagai et al. 2007 Arnaud et al.|2010[ |, 
and <5 = 0.12 is a modification which better captures the vari- 
ation of mass scaling with radius in the REXCESS sample 
( [Arnaud et al.|2010[[Bohringer et al.|2007^ . 



file ( [Navarro et al.|1997| [Nagai et al.|2007[ ). We discretize 
this profile by a volume-weighed average over our radial 
bins, which makes for a more meaningful comparison with 
the discretized universal pressure profile we have defined, 
i.e. {Pfe}. Our measurements are in good agreement. 

Table |2] shows the level of confidence for our vari- 
ous pressure measurements. The difference between the 
standard and modified scalings of Pc with mass is very 
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Table 3. The normalized covariance matrix of the universal pressure profile for all MCXC clusters (blue or top-right), as well 
as only the resolved ones (red or bottom-left). The modified self-similar model is used for these measurements (i.e. S = 0.12 in 
Equation dQl). To construct these matrices, let Cp" be the covariance matrix for analysis done on all MCXC clusters. Construct a 

diagonal matrix A such that A| j, = \/[Cp ]k,k, where k runs over the different radial bins. We define the normalized covariance 
matrix via D''" = [ A^"] ^ ^ Cp [A^"] ^ ^ , which normalizes all diagonal elements of Cp to one. By the same construction, let D'^'^'' be 
the resulting normalized covariance matrix for analysis done on resolved MCXC clusters. The blue or top-right (red or bottom-left) 
numbers in this table denote the off-diagonal elements of D^" (D''""), respectively. 



blue 



Bin 1 Bin 2 Bin 3 Bin 4 Bin 5 



Bin I 



Bin 7 Bin 8 



Bin 1 


1.000 


-0.727 


0.282 


-0.082 


0.029 


-0.022 


0.011 


-0.001 


Bin 2 


-0.612 


1.000 


-0.671 


0.214 


-0.051 


0.005 


-0.003 


0.006 


Bin 3 


0.277 


-0.497 


1.000 


-0.633 


0.144 


0.000 


-0.023 


0.018 


Bin 4 


-0.069 


0.232 


-0.445 


1.000 


-0.598 


0.109 


0.009 


-0.014 


Bin 5 


0.033 


-0.039 


0.147 


-0.434 


1.000 


-0.594 


0.117 


-0.004 


Bin 6 


-0.024 


-0.007 


-0.013 


0.091 


-0.451 


1.000 


-0.597 


0.156 


Bin 7 


-0.001 


-0.023 


-0.041 


-0.015 


0.108 


-0.442 


1.000 


-0.697 


Bin 8 


-0.011 


-0.043 


-0.044 


-0.075 


-0.053 


0.085 


-0.476 


1.000 



Table 2. Level of detection for various pressure measure- 
ments. (5 = (0.12) corresponds to measurements presented 
in Fig. |3a| l |3b[ |, respectively. The null chi-squared is given by 
Xo = P Cp P, where P are the best fit pressure measure- 
ments and Cp is their associated covariance matrix. The level 
of detection is calculated for 8 degrees of freedom, i.e. number 
of radial bins. 



Measurement 



xl 



Detection (cr) 



All MCXC clusters, 5 = 
All MCXC clusters, 5 = 0.12 
Resolved MCXC clusters, <5 = 
Resolved MCXC clusters, <5 = 





259.3 


15.1 




262.2 


15.2 





115.6 


9.5 


0.12 


118.6 


9.6 



small. However, the significance of detection reduces from 
15.1(7 to 9.5(7 if we limit our sample to the resolved clus- 
ters. This may seem surprising because, looking at Fig. 
[3] the uncertainties are similar in both cases and the best 
fit pressure values are even consistently higher in the case 
of resolved clusters. Note, however, that Fig. [S] does not 
compare the off-diagonal elements of the covariance ma- 
trices — i.e. correlation between different bins. In fact, 
the extra statistical information coming from unresolved 
clusters is encoded almost entirely m the off-diagonal cor- 
relations. We refer the reader to Appendix [E| for the full 
covariance matrix and numerical values of {Pfc}. To get 
some sense for the nature of correlations, however, we 
have shown a normalized version of Cp for both samples 
in Table [3] Nearby bins are anti-correlated in both cases, 
but more so for the sample containing all clusters. The 
extra information contained in these anti-correlations can 
be quantitatively described by examining the eigenvalues 
and eigenvectors of the covariance matrix Cp. We refer 
the reader to AppendixlDlfor a detailed discussion of this 
point and state the results here. In the case of the resolved 
clusters, the eigenmodes with the three largest eigenval- 
ues are responsible for most of the contribution to Xo- 
For the whole MCXC sample, however, all eigenmodes 
contribute more or less equally. Moreover, eigenvectors 
corresponding to larger eigenvalues carry most of their 



weight from the inner bins. Therefore, in the case of re- 
solved clusters, mostly the inner bins are contributing to 
the signal, whereas for the whole MCXC sample, there 
is also contribution from outer bins. This analysis reas- 
sures us that even unresolved clusters contribute to the 
tSZ signal in the outskirts of the ICM. 

Although the unresolved clusters add to the tSZ sig- 
nal in the outer bins, one expects numerical uncertain- 
ties associated with them. This is especially worrisome 
for those on sub-pixel scales, where certain approxima- 
tions, such as a spherically averaged pixel transfer func- 
tion, break down. In order to get an estimate for how 
large such uncertainties are, we performed our analysis 
on all MCXC clusters using higher resolution WMAP sky 
maps (Nside ~ 1024). The result is shown in Appendix 
B2| For all radial bins, this discrepancy is at most at the 
1(7 level and is random in nature. 

Our pressure measurements are also affected by the 
uncertainty present in determining masses of clusters. In 
Appendix [Cl we have investigated this issue by consider- 
ing 62 MCXC clusters which are also in the Early Release 



SZ (ESZ) sample ( jPlanck Collaboration VIII|2011b| ). The 
ESZ mass estimates are systematically higher on average 
by about 12 percent. This results in systematically lower 
pressure measurements (where there is actual signal), but 
it is only at the la level (see Appendix [c] for details). We 
will return to this issue in Section [4.2[ because this effect 
is no longer small when determining gas mass fraction. 

The results of our analysis as applied to the cluster 
subsamples introduced in Table IT] are included in Ap- 
pendix [El Because of the large statistical uncertainties, 
comparing the pressure profile of different mass bins is 
not terribly illuminating. We will, however, discuss the 
implications for gas mass fraction in the next Section. 

4.2 Gas Mass Fraction 

The density of gas pg (r) in the a**^ cluster with temper- 
ature profile T^'^'lr) takes the form 



prir) 



^enipPr (r) 



(27) 



kBT(°^'>{r) ' 
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where nip is the proton mass and fie — ^^q-j- ~ 1.14 is 
the mean molecular weight per free electron for a cosmic 
hydrogen abundance oi X ~ 0.76(_] As it has been the 
case for the electron pressure profile (see Equation (|6|), 
we assume a universal temperature profile 

T<'''(r) = ri"'T(r/r<;"'), 

which in turn implies 



(28) 



» 



(r) 



U)nr/rn 



T(r/r, 
(a) _ fJ-empP, 



{<i)\' 



(a) 



knT, 



(a) 



(29a) 



(29b) 



The volume-averaged gas density at radius r takes the 
form 



pT'i<r) 



i-Kr'^ jQ 



3pi''\r/ri'^'>y' 



pg (r )4:Tvr dr 



T(x 



X) 2 , 

X ax. 



(30) 



Given that we have considered radial bins throughout 
which P is constant, we shall also bin the temperature 
profile: 



Tfc 



1 



fK-^Li) 



T(a;)a; dx, 



(31) 



where n^ is the value of thefc"^ radial bin in units of r^ , 
with no = (see Section |2.2|.p| It then follows that 



4"'(<0 = pi"'EV',(r/r 



(«)\ 



fc, 



(32) 



where Nh is the total number of radial bins and 
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radius r is equal torj 



m(4"V/i?i"0,^w^ 



Mi") (< r) = •"^-^;/;;^ ' m^' , (35a) 



m{x) = ln(l + a;) ■ 



1 + x' 



(35b) 



where R^ and A/^ were defined in Equation ([8|. Also, 
it follows from Equation ([sl that 






35(") ■ 



(36) 



Therefore, knowing c^ (for any A) determines 5'"\ or 
equivalently pi" . We estimate the concentration param- 
eter from the C200 — M200 relation of |Bhattacharya et al.| 
( [20TT| : 

(37a) 






+ 0.53 , 
(37b) 



where D(z) is the linear growth factor normalized to 1 
at 2 = 0. fj As was the case with the temperature of 
baryons, we similarly bin pm (r) 

W _Mi°>«r<"')-Mi°'«rt\) 



Pm,k ^ 



f {[-i"']^-K-\]^} 



(38) 



where as before r 



(a) 



averaged matter density pin up to radius r is 



rikr. 



(a) 



Finally, the volume- 



p!:'«0=EVl.(r/r("')p: 



m.k ' 



(39) 



'^[nl~nl_,] iik<k. 



Vk{x) 



1 



[^'"<] 



.0 



if fc = fc. -fl 
if fc > fc, + 1. 



(33) 



Here k, is an integer defined via Uk, < x < nk,+i. We 
assume the total matter density of the a**^ cluster to be 
of the NEW form: 



J«) 



P\^>{r)^ 



■/ri"'(l + r/r<"))^ 



(34) 



Defining 4"' = E^"'/^^"' and ^t"' = pi^Vpc^^^-), it 
may be checked that the total mass enclosed within a 



where {k, being defined as above) 



^ H - ^l-i] 







x'-nl] 



if fc < fe. 

if fc = fc* -^- 1 

if fc> fc. + 1. 



(40) 



The average gas mass-fraction up to radius r in the a**^ 
cluster takes the form: 



n:u<r) 



P^'"«r) 

n [ 



p)^>(<r) 

(a) N^ 



(41) 



* The free electron number density isrie = nff + 2n[je, where 
riH and nue are the Hydrogen and Helium number density. 
The cosmic hydrogen abundance is X = n[j/{n[j + Anff^)- It 
then follows that pi, = mplriH+AnHe) = v- .'^ rte = y+i k °T ' 
where rUp is the proton mass. 

^ We will discretize all continuous profiles over our radial 
bins, because our pressure measurements are discrete by con- 
struction. In the case of temperature, it might seem more nat- 
ural from Equation l |30[ l to discretize 1/T instead of T. We 
have checked that the difference between these discretization 
schemes is insignificant. 



10 



Here we have adopted the notation of|Bhattacharya et al.| 



(20111 



TT 



DiM 



A good approximation to D{z) is given by D{z) — „ ,„. , 
where (see [Carroll et al.|1992| 

6fl,n(z) 

Di{z) ~ — ^x 

2(1 + z) 

[n,n{z)'^^^ - ^Aiz) + [1 + n,n{z)/2] [1 + nA{z)/70]y . 

Here Cm(z) = Qrnil + z)^/E(zf and J^aCz) = ^a/E{z)^. 
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In order to make a meaningful comparison with the uni- 
versal gas mass-fraction, we average this quantity over all 
clusters 



(42a 



fgas{< X) -= 


= — V f*" 
a — l 


« xr'f 


') 






= A{x)^V, 
fc=i 


(x)Pfe, 






A{x), 


" a = l Pm 


pi-' 
(< xr'^ 


~) 



(42b) 



Since fgas{< x) is a linear combination of {Pfc}, it is a 
Gaussian random variable with mean and variance 



{fgasiK X)) = A{x)Y,Vk(x)fk 
fe = l 

^}ga.s«^) = {[fs'^si< X) - {fgas{< x))f) 

^Aixf Yl V,.{x)Vk>{x)[Crh 



(43) 



(44) 



The averaging scheme we have adopted in Equation ( 42 1 



may seem arbitrary. One could, for instance, assign dif- 
ferent weights to different clusters. If Pc and Tc scale 
similarly with mass, different averaging schemes would 
differ by a negligible amount. This is because the only 
variation in fgas amongst different clusters would be due 
to the scaling of csoo with mass, which is fairly mild. As 
a result, given the temperature profile we have adopted 



(see Equation (45 1), we use the standard self-similar scal- 
ing of Pc with mass (5 = in Equation ([9| below) to 

compute fgas- 

We use the average temperature profile of |Vikhlinin| 
[eFahlpOOet : 

,(a;/0.045)^-^+0.45 1 



T{x) = 1.35- 



rp(a) 

5 keV 



(j:/0.045)i-9 + 1 [1 + (3;/0.6)2]'''« ' 

(45a) 

1/1.51 

(45b) 



MtolEiz) 



3.41 X 10i4/!.-iM, 



This is an approximation to the averaged profile of about 
a dozen low-redshift X-ray clusters, with measurements 
obtained for r < R500 ( Vikhlinin et al.|2006 1. The scatter 
about the mean profile is about 15 percent. The assump- 



tion of universality (i.e. Equation (28l) may be easily re- 
laxed if temperature measurements for individual clusters 
are available. In the case of our present work, however, 
this option is not viable since a large cluster sample is 
required to compensate for WMAP's insufficient sensi- 
tivity. To get an estimate for how this assumption affects 
our fgas measurements, consider an average lOOe percent 
scatter about the universal profile. For a relatively large 
sample of clusters, as is the case with our measurements, 
changing the individual temperature profiles by ~ lOOe 
percent would introduce a systematic uncertainty of or- 
der e^ in /gas -F^ Given the large statistical uncertainties 
in our pressure measurements, this effect is small. 

12 This is because fgas oc Y,a 1/7"^"', where T'") is the tern- 
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(a) Gas mass fraction for all MCXC clusters. 



— Mean Profile of Resolved MCXC Clusters 

— - Systematic Uncertainties 

— Cosmic Baryon Fraction 



(b) Gas mass fraction for resolved MCXC clusters. 

Figure 4. Gas mass fraction of all MCXC, as well as the re- 
solved subsample. A cluster is considered resolved if its first 
radial bin subtends a solid angle larger than the effective beam 
area of the W frequency channel (see Section 3.2). The stan- 
dard self-similar model of pressure is assumed for both mea- 
surements (5 = in Equation (JOI). The black solid curves 
show the average gas mass fraction, computed using Equa- 
tion Bsl. The shaded areas represent the standard deviation 
in the measurement of fgas as given by Equation ( |44[ |. The 
dashed black lines show the expected systematic uncertainty 
about the mean gas mass profile, mostly due to cluster mass 
estimates. The cosmic gas mass fraction is obtained by fit- 
ting ACDM to WMAP 9-I-SPT+ACT data a nd is equal to 
ftt/nm = 0.165 ± 0.005 pinshaw et al.|2012| l. 



Systematic uncertainties associated with estimating 
masses of clusters have a more significant effect on gas 
mass fraction. As mentioned in the previous section, we 



perature of the a cluster. Let {ea} be realizations of a gaus- 
sian random variable with zero mean and standard deviation 
e. Changing the temperature T'"' — ^ T''''(1 + ea) is equiva- 
lent to J2a l/'T'"' -^ Ea l/r*"' (1 - ea -l-e^ -I- • • ■ )• When T(") 
doesn't change drastically from cluster to cluster, the contri- 
bution of the term linear in ea is on the order of yNe, while the 
second order term contributes about Nt'^. When A'^ > 1/e^, 
the e^ term dominates. 
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have investigated this issue in Appendix [C] by consider- 
ing 62 MCXC clusters which are also in the ESZ sample 
(iPlanck Collaboration Villi 1201 Ibl). The ESZ mass es- 



timates are systematically higher by about 12 percent, 
which causes lower pressure measurements at the la 
level. Because of the scaling of temperature with mass 



(Equation (45 1), this decreases fgas by about 20 percent 
(see Fig. |C1[ ). Repeating our analysis with Planck CMB 
data is expected to reduce statistical errors significantly 
(see Section p] and Planck Collaboration V (20131). In 



that case, fgas measurements would be solely dominated 
by systematic uncertainties associated with determining 
masses of clusters. 

Fig. H] shows the result of our analysis, as applied 
to all MCXC clusters, as well as just the resolved ones. 
Given the large statistical and systematic uncertainties, 
we extrapolated the temperature profile of [Vikhlinin] 
et al. (20061 out to R2oo- The black solid curves show 



the average gas mass fraction, computed using Equation 
( 43 1 . The shaded areas represent the standard deviation 
in the measurement of fgaa as given by Equation (44 1. 



The dashed black lines show the expected systematic un- 
certainty about the mean gas mass profile (~ 20 per- 
cent), mostly due to cluster mass estimates. Considering 
both statistical and systematic errors, our results are fully 
consistent with the cosmic baryonic fraction up to -R200- 
Given the large error bars, accounting for all baryons in 
stars does not change this conclusion. 

Fig. [5] shows the dependence of gas mass fraction 
(up to -R500) on the cluster subsamples which we have 
binned according to mass (see Table [I|. A general trend 
of increasing fgas with M500 can be seen. Due to the 
large error bars, however, we cannot make any statisti- 
cally significant statement about this dependence. Fol- 



lowing Planelles et al. (20131, we compare our measure- 



ments with two different observational samples: a com- 
bined sample of 41 clusters and groups from |Vikhlinin| 



et al. (20061, Arnaud et al. (20071 and Sun et al. (20091 



(V06-I-APP07-I-S09), shown as the green region, and the 
sample obtained from the combination of the data by 



Zhang et al. (20111 and Sun et al. (20091 (Z11-HS09), 



shown as the yellow area (see Table 1 of Planelles et al. 
(2013l). The black stars show fgas obtained from hy- 
drodynamical simulations carried out by [Planelles et aL] 
p013 ), which include radiative cooling, star formation 
and feedback from supernovae and active galactic nuclei. 
Considering both statistical and systematic errors, our 
measurements are consistent with both the observational 
and numerical results. 



5 DISCUSSION AND FUTURE WORK 

In a companion paper, we will repeat our entire analysis 
using the recently released Planck CMB data. To get an 
idea for how much our measurements will improve, we es- 
timate here the pressure covariance matrix Cp expected 
from Planck. We consider the six Planck-'HFl channels, 
which have central frequencies 100, 143, 217, 353, 545, and 
857 GHz, at HEALPix resolution Nside = 2048. We as- 
sume a homogeneous detector noise constructed by aver- 
aging the noise variance of all pixels for a given frequency 
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(a) Binning the entire MCXC sample according to mass, as 
shown in Table llal 
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(b) Binning the resolved MCXC sample according to mass, 
as shown in Table llbl 

Figure 5. Dependence of fgas on the mass of clusters. All 
mass bins are defined in Table [T] The black dots show fgas 
measurements up to -R50O1 with black vertical bars denoting 
statistical errors. The red vertical bars show systematic un- 
certainties expected due to cluster mass estimates. Following 
[Planelles et al. ] | [2013[ l, we compare our measurements with 
two different observational samples: a combined sample of 41 
clusters and groups fro m[Vikhlinin et al.[ pOOM , [Arnaud et al.[ 
( [2007| l and |Sun et aLJ ( [2009| l (V06+APP07-I-S09), shown as 
the green region, and the sample obtained from the combina- 
tion of the data by [Zhang et al.[ ( [201lt and [Sun et"al] |[2009[| 
(Z11-I-S09), shown as the yellow area (see Table 1 of |Planelles[ 
[et al.[ (|2013i)). The black stars show fgas obtained from hydro- 
dynamical simulations carried out by [Planelles et al.[ poTsl , 
which include radiative cooling, star formation and feedback 
from supernovae and active galactic nuclei. 



channel: 



[Cn] 



2 r r 



(46) 



where nioo = 50.6 /xK, 71143 = 20.1 /xK, 71217 ~ 27.1 /iK, 
71353 =0.1 mK, 71545 = 28.1 mK, 71357 = 27.9 mK. We 
assume Gaussian instrumental beams with FuU-Width- 
Half-Maximum of 9.5,7.1,5.0,5.0,5.0,5.0 arcmins, for 
the six HFI channels in order of increasing frequency 



(Planck HFI Core Team VI 20111. Finally, we assume 
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no masking (i.e. M = 1) but remove all clusters that are 
masked out from our templates. 

Because we have assumed a homogeneous detector 
noise and no masking, the matrix D introduced in Equa- 



tion ( 23 1 can be inverted analytically, resulting in the 
pressure covariance matrix: 



[Cj 



E E E(4" 



f(-)\ 






ImC,^ 



('•-) 



Ci 



1 + Ce^t 



4i'^)\, cM 



E(*t')^".e 



M\. cM 



E(4'')'"^e; 



(47) 
L 1^ J L i' 

where 

(48a) 

(48b) 

(48c) 
Here Apix denotes the pixel area and all other quantities 





{Bt.Wi 


iVf' 


^i 


-Ed^' 


v(^' 


— -^pixf^u ■ 



are defined in Sections 12.31 and 12.41 

We generate our templates t^ using the resolved 
MCXC clusters which are not masked (total of 122 clus- 
ters), with a standard self-similar pressure-dependent 
scaling (5 = in Equation ([9|). In order to make sure 
this estimate is reasonable, we computed the same quan- 
tity with WMAP9 data, using the Q, V, and W channels. 
Assuming the best-fit pressure values P remain the same, 
this leads to a null chi-squared of Xo = 124.551, which is 
reasonably close to the actual value Xo = 115.626. (We 
use the same 8 radial bins as for our pressure measure- 
ments, i.e. there are 8 degrees of freedom here.) Esti- 
mating the covariance matrix for Planck using the same 
best-fit pressure values, we obtain Xo ~ 66154.8. There- 
fore, assuming that the signal does not change, we ex- 
pect the statistical uncertainties to reduce by a factor of 
~ ^66154.8/124. 551 = 23.9. This is a significant im- 
provement, which will allow us to consider finer bins and 
possibly probe the ICM pressure to larger radii. Fig. |6] 
compares the expected error for different bins with those 
of WMAP9. 

We have not addressed the issue of point source con- 
tamination so far. In our framework, it is not feasible 
to fit locally for contribution of point sources, given the 
large number of clusters. We did try to account for them 
by assigning a constant absolute luminosity per frequency 
channels to all clusters. The results, however, change only 
by a negligible amount. We will consider a more detailed 
modelling of point source contamination in future work. 



6 CONCLUSIONS 

We have introduced a statistically-optimal and model- 
independent framework for extracting the universal pres- 
sure profile of the hot gas in the intracluster medium. 
The thermal Sunyaev-Zeldovich effect makes this possible 
because it is linearly proportional to the integral of the 
electron pressure along the line of sight. We use the prin- 
ciple of maximum likelihood to find best-estimate values 



WMAP9 Noise 

Estimated Plancl< Noise 

Battaglia et al. (2012) 

Expected Signal: Nagai et al. (2007) 




Figure 6. Comparison of the statistical noise expected from 
analysis with Planck CMB data vs. WMAP9. The blue data 
points are the estimated noise expected from repeating our 
analysis with Planck CMB data (see Equation EtI). The red 
points correspond to the same quantity for WMAP9 measure- 
ments, applied to the resolved MCXC clusters with standard 
self-similar pressure-dependent scaling (<5 = in Equation 
^n). The green and black curves show the expected pressure 
signal from simulated clusters of [Battaglia et al.] | |2012| and 
[Nagai et al.|||20()7ll, respectively. 



for the radially binned values of the pressure profile, as 
well as the full covariance matrix governing their uncer- 
tainties. Once reformulated in the proper mathematical 
framework, the main technical challenge is solving a very 
large system of linear equations, which we do numerically 
by employing the conjugate gradient method. 

We applied our methodology to WMAP9 data and 
various subsamples of the MCXC catalogue. In the case 
of all MCXC clusters, we extract the pressure profile with 
a high accuracy at ~ 15(t confidence, with possible sys- 
tenratics uncertainties dragging our detection down to 
~ 14ct. We also considered a subsample of the MCXC 
clusters completely resolved by the W frequency channel 
of WMAP, resulting in a ~ 9a detection. In an upcoming 
companion to this paper, we apply the same methodol- 
ogy to the recently released Planck CMB maps. An esti- 
mation of the pressure covariance matrix expected from 
Planck suggests that the current signal-to-noise will im- 
prove by a factor of ~ 24. 

Assuming a temperature profile motivated by X-ray 
observations, we computed the average gas mass fraction 
as a function of radius. We argued that systematic uncer- 
tainties associated with estimating mass of clusters could 
have a drastic effect (~ 20 percent) on gas mass frac- 
tion. Considering both statistical and systematic errors, 
our results are fully consistent with the cosmic baryonic 
fraction and the expected gas mass fraction in halos, up 

to 7?200- 

We also made a first attempt at studying the depen- 
dence of gas mass fraction on the mass of clusters. Due 
to the large error bars, we cannot make any statistically 
significant statements about this dependence. Neverthe- 
less, our measurements suggest that gas mass fraction 
increases with the mass of clusters, which is consistent 
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with findings from X-ray measurements (Viklilinin at al. 
20061 lArnaud et al.|[2057l |Sun et al.||2009[ [Zhang et al. 



20111 and numerical simulations (Planelles et al.|2013l 
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APPENDIX A: TECHNICAL DETAILS 
Al Likelihood Function 

Given an underlying temperature field rj 



5T(n) ^YIY. STimYi^in), 



the value of the discretized temperature map at pixel i and frequency band v is given by 

Imax t 



"^^^ =Y. Y. STe^Be.WiYe^in,), 



(Al) 



(A2) 



where B^^ is the isotropicized beam transfer function for the mode I and frequency channel u, and We is the isotropicized 
pixel transfer function. In case of the primordial anisotropies, because STf^, are linear functionals of ST^ (n), they are 
correlated Gaussian random variables with zero mean and two-point function 



[Cs]i^,i>i,' = {5Ti„ 5Ti,^,) = 2J 



2£+l 
47r 



CiBt,Be^,WtPi(h,-h,,), 



(A3) 



where Pi is the ^-th degree Legendre polynomial and we have used {ST^^ 5T^^,) — CiSuiSmm'- 

Finally, the log-likelihood probability — x^/2 of jointly measuring the CMB temperature values {^Ti^j^g^" given 

the tSZ contribution {ST^J^}'^^^" and the primary CMB fluctuations {STfJ^Yi^^" is 

- lx\{ST..}:ii:\{STi^}-ii:,{STF^}-i^;) = 

- ^<5tP^^Cs"'<5T''* - hsT - ST''-^ - ST'^^fC^^^ST - ST'''^ - ST^^), (A4) 

where Cn is the noise covariance matrix and Lp {Lv) is the set containing pixels (frequency channels) we wish to use 
in our analysis, rj After integrating over all possible primary fluctuations ST^ , which can be done analytically, the 
log-likelihood takes the form 



xiiST^.V^itliSTrfV^it) = - =,iST - ST-y C-\ST - ST^^) 



(A6) 



with 



This proves Equation ( 10 1 



C = Cs + C, 



(A7) 



A2 Masking 



All quantities used (and not defined) here are introduced in Sections 2.3 and 2.4 



Let Lp (Lp respectively) denotes the set of all masked (unmasked respectively) pixels, so that L = LpU Lp contains 
all pixels on the sky. In what will follow, Csf (Cs) will denote the signal covariance matrix defined on L (Lp). The 
same notation will be used for the noise covariance matrix. Consider now the full covariance matrix Cf — Csf + CNf . 
It is related to C = Cs + Cn via C — Cp"^CfCp, where [Cp]ij^^jv' = PijS^^i is a projection matrix with i £ L and 
j & Lp, such that all components of P are zero except Pa = 1 for i £ Lp. Construct from Cf another matrix Cf with 

— X for all i £ Lp. Then, it can be shown from the definition of the inverse of a 



the same entries, except that [Cf] 



^* We follow the HEALPix conventions for spherical harmonic transforms I Gorski et al. 12010 i 
^^ The matrix notation used here is explicitly: 

where we are considering Cg as a matrix whose rows and cokimns are labeled by iu. 



<5TPA^Cs-15TPA 



(A5) 
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matrix that! 



C"^ = lim Cp^Cf Cp. 



(A9) 



Because [CNf]i,/,i'iy' = [Nv]ii'S,^^i , we may construct Cf as follows: [Cf]ii/,i'^' = [Ni/Jii'^^i/'+Csf , where N^ has the same 
elements as N,^ except [N^Jw = x for all i G Lp. Since [N^Jii' = n^j^Sa' is diagonal, it follows that [N^ ]iii — l/[Ny]ii', 
which then implies 



lim N„ ^ = mn;:^m, 



where M = PP"^ is the masking matrix defined in Section 2.4 Moreover, we show in Appendix A3 thatrj 



[^f liiy.i' ii' 

Using the above two equations we find 



n:'5.,., - n:'b.s^/^ 1 + s^/^ y b,n:'b, s^/^ s^/^b.,n;/ 



J2 b„n; b„ ) s^ 



lim [Cf liuyv' — [G^^^'liii 



(AlO) 



(All) 



(A12) 



where G^_^/ is given by ( |23b[ ). 

Finally, let [V/Ji,/ and [W/]i,y be two vectors defined on every pixel on the sky (i.e. i £ L). Also, let V and W be 
the corresponding vectors defined only on the unmasked pixels: V = Cp"^V/, W = Cp"^W/. Then 

(A13) 
(A14) 

(A15) 



V'^C'^W = lim V>^CpCp^C7^CpCp^W/ 

a;— S-oo 






i,i' GL h',!^' ^L,^ 



where we have used the fact that [CpCp-^Ji^^/^/ = MaiSy^i, as well as M^ = M, which combined with ( |A12[ ) implies 
MGt,,^/M = G^^v'- This justifies the equality between the corresponding equations in (151 and (22 1. 



A3 Covariance Matrix Re-loaded 

All quantities used (and not defined) here are introduced in Sections |2.3| and |2.4[ Also, all matrices are defined on the 
entire sky. 

Using the definition of matrices B^ and S as given in Section [2. 4| it may be checked that 



where we have used 



It then follows that 



Api:cy^Ylm{hi)Yl,jn'{ni) ^ / Yi™,(n)y'i'm'(n)rf^W = 5;;'5„ 



iGL 



Let us now prove the following 






where 



G,y = ^Z^5,,y - CpN;'B.Si''''D-'s'''"B,yN,:/, D = 1 + S^''' i ^ B.N^^B, j S'^^ 



(A16) 
(A17) 

(A18) 
(A19) 

(A20) 



16 



As a simple example, consider the case where there are only two pixels, one frequency channel, and one of the pix- 



els is masked out: Lp = {1}, Lp = {2}. In this case Cp = P 
-1 1 



, C = [Cfhi, Cf 



[Cf]il [Cf]i2 
[Cf]21 X 



Cf 



^[Cf]ll-[Cflf2 



X -[Cf]l2 

-[Cf]l2 [Cf]ll 



lim Cp^Cf ^Cp = lim 



Equality jA9| can now be easily verified: 
X 1 



:C-l. 



II2 L'-fJli 



and 



(A8) 



^^ Note that the proof presented in Appendix A3 proceeds in the exact same fashion when N,y is replaced by N^. 
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That this is true may be easily checked: 

= [B,.SB,.N;i].fe + [N,.N;i].fc5,p - [B,S^/'D-^S^/'CfN;i].fe 
B.S^/^l + S^/^ ( E B.N^^B.'j S^/4 D-^S^/^B.N;1 

+ [B,S'/^D-'s'/^B,N;'],fc 

= [B,SB,.N;i],fe + S.kS^p - [B,S'/'DD-'s'/'B,N;i],fc 

where we have used S'^'^^S''"^^ = S, which may be checked using (A17|. 



(A21) 



A4 Fitting the Monopole and the Dipole 

Let us briefly discuss how any possible residual monopole and dipole CMB components can be accounted for in our 
framework. We denote monopole and dipole contributions by ^T^^" and 5T^^^ , respectively. They take the form 



5T'-\h) = aooyoo(n), 



5T'-^{h) = awYia{h) + 2Re(aii)Re(Yii(n)) - 2Im(aii)Im(yii(n)). 



(A22a) 
(A22b) 



These components should be added to the tSZ signal: 5T^'^{h,v) -^ 5T™{h,v) + ST'^-°{h) + ST'^-^{n). This may be 
conveniently done by making the following definitions: 



iVb+i = ffloo, Jr]v^+2 = aw, 



^+3 = Re(aii), Pjv^+4 = Im(aii) 



(A23a) 



4,+i(^) = >'oo(n), <V2(^) = yio(n), 4;V3W = 2Re(yii(n)), 4;V(^) = -2Im(ni(n)), (A23b) 

where Nb is the total number of radial bins. The statistical machinery developed in Section [23] now goes through exactly 



the same way, except that Nb ^ Ni, + 4. Once the matrix a is found (see (15l), which is now {Nb + 4) x {Nb + 4) 



dimensional, the pressure covariance matrix becomes the restriction of its inverse to the bins of physical interest: 

[Cp]fcfc. = [cx-\k' where k,k' &{1,..., Nt}. 



APPENDIX B: ROBUSTNESS TESTS 

As it was shown in Section [2. 4| the most important part of our analysis is solving a linear equation of the form Ace = b, 
where A is a very large (~ 10" x 10^) matrix. In order to be certain that our numerical methods are correct, we perform 
two tests. 



Bl Simulations 

We create sky maps with known tSZ amplitudes (i.e. the quantities of interest Pfe and other parameters such as 
the monopole and dipole anisotropics) and see if the outcome of the pipeline matches with what is inputted. More 
specifically, we generate A'' random realizations of the CMB primary anisotropics, add the tSZ signal with known 
amplitudes Pfc, and finally add random detector noise. The outcome of every fitting procedure is the set of values Pj, , 
with i £ {1, . . . ,N} denoting the i-th simulation, and the covariance matrix [Cpjfej;/ (see Section 2.3 1. The covariance 



matrix [Cp]^^^/ does not change from one simulation to the other since it only depends on the tSZ templates, and the 
detector noise and primary CMB covariance matrices. If Pj, really are realizations of a gaussian random variable with 
mean Pj, and variance [Cp]fe,fc, then their mean P^ — ;^ X^j^i Pj, should converge to P^ as A'^ becomes large. More 

specifically, the expected error in determining the true value of the mean is 



the best estimator of the variance a^ = 

V< {^l - [Cp],,fe)2 > = y^[Cp],,fe. Fig. 



< (Pfc - Pfe)2 > = ^[Cr]k,k/N. Similarly 
X]i=i(lPfe ~ Pfc)^ should converge to [Cp]^,^', with an expected error of 
31 shows the results of our simulations for a few templates (i.e. values of 
k). As it can be seen, all estimators converge to the values computed by our pipeline. 
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Figure Bl. Testing the pipeline with random primary CMB+noise simulations. The shaded area represents the expected error 
in the quantity of interest (see text for more details.) 

B2 Different Resolutions 

We performed our analysis on all MCXC clusters using WMAP7 sky maps at two different HEALPix resolutions of 
9 {Nside ~ 512) and 10 {Nside = 1024). The modified self-similar scaling of pressure with mass is used for these 
measurements (i.e. 5 = 0.12 in Equation Q). Fig. |B2| shows the results. The null chi-squared for the Ngidf, = 512 
(Nside = 1024) measurement is Xo ~ 246.8 (xo = 272.1), corresponding to a 14.68a (15.50(t) detection. 



APPENDIX C: EFFECT OF UNCERTAINTY IN MASS OF CLUSTERS 

Our entire analysis depends crucially on the self-similarity length/mass scales of clusters. Therefore, it is important to 
investigate how our measurements are affected by the uncertainty present in determining masses of clusters. Planck| 



Collaboration V (20131 use 62 clusters from the Early Release SZ (ESZ) sample (Planck Collaboration VIII 



2011b I 
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Figure B2. Results of our analysis of WMAP7 sky maps at two different resolutions. 



which also belong to the MCXC catalogue (Plan ck Collaboration XI|2011 \. To get an idea for the degree of uncertainty 
present in mass estimates, we compare ESZ and MCXC masses of these clusters. This is shown in Fig. |Cla[ The ESZ 
mass estimates are systematically higher on average by about 12 percent. To investigate how such systematics affect 
our pressure measurements, we randomly changed masses of all MCXC clusters according to the distribution in Fig. 
Cla (We used the standard self-similar scaling to create our templates, i.e. we set 5 = in Equation (|9|.) The resulting 
pressure profile is shown in Fig. ( |Clb[ ). In the first three bins, where there is signal, the pressure values decrease 
systematically. This difference, however, is at most at the la level. As shown in Fig. |Clc[ this is no longer the case for 
gas mass fraction, which decreases by about 20 percent on average. 



APPENDIX D: ALL VS. RESOLVED MCXC CLUSTERS 

Let us provide a quantitative explanation for how statistical information is lost when unresolved clusters are not 
accounted for. The discussion that will follow is based on measurements presented in Fig. |3b| i.e. the modified scaling of 
pressure with mass {5 = 0.12 in Equation Q). The same analysis for the standard self-similar scaling (i.e. measurements 
presented in Fig. 3b I gives the same results. Let A„ and T„ denote the eigenvalues and eigenvectors of Cp, respectively. 
We choose the labels n such that Ai < A2 < • ■ • < AjVj^, where Ai';, = 8 is the total number of radial bins. Since Cp 
is a positive symmetric matrix, its eigenvalues are positive and its eigenvectors are real. The null chi-squared can be 
re-expressed as 



2 
Xo 



E( 



i^T„)VA„ 



(Dl) 



Fig. 



Dl 



shows the eigenvalues A„ and the contribution (P"^Tn)^/A„ of the different eigenmodes to xo- In the case of 
the resolved clusters, the modes with the three largest eigenvalues are responsible for most of the contribution to xo- 
For the whole MCXC sample, however, all eigenmodes contribute more or less equally. Eigenvectors corresponding to 
larger eigenvalues c arry most of their weight from the inner bins. To see this, we have plotted the components of all 
eigenvectors in Fig. D2 We denote the fc**" component of the eigenvector Tn by T^^k- As in the text, k — I . . . Nt, labels 
the radial bins around clusters with fc = 1 and k = Nb corresponding to the inner and outer-most bins, respectively. 
It is clear from Fig. |D2| that for eigenvectors with larger eigenvalues, the components corresponding to the inner bins 
dominate, and vice versa. Therefore, in the case of resolved MCXC clusters, the fact that most of the contribution to Xo 
comes from modes with the three largest eigenvalues indicates that mostly the inner bins are contributing to the signal. 
Whereas for the whole MCXC sample, there is also contribution from outer bins. This analysis reassures our intuition 
that even unresolved clusters contribute to the tSZ signal in the outskirts of the ICM. 
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(a) ESZ vs. MCXC mass estimates of 62 clusters common to 
both catalogues (Planck Collaboration XI|2011 1. 
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(b) Effect of cluster mass uncertainties on the universal 
pressure profile. 
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(c) Effect of cluster mass uncertainties on gas mass profile. 

Figure CI. Effect of systematic uncertainties associated with the mass of clusters, (a) Difference between ESZ and MCXC mass 
estimates of 62 clusters common to both catalogues ( [Planck Collaboration XI|2011[ |. (b) Effect of changing MCXC masses on the 
universal pressure profile. The entire MCXC sample is used and the standard self-similar model is assumed (5 = in Equation 
(J9|). The blue data points show the result of our measurements using MCXC mass estimates. The red points show measurements 
for which the MCXC masses are randomly changed according to the distribution shown in (a), (c) Gas mass fraction, computed 
using Equation ( |43[ l, corresponding to pressure measurements in (b). The shaded areas represent the standard deviation in the 
measurement of fgaa as given by Equation (J44l. 



APPENDIX E: PRESSURE MEASUREMENTS: THE EXACT NUMBERS 

Here we report exact numbers corresponding to the measurement of P, i.e. the best fit universal pressure value P and its 
associated covariance matrix Cp. As before, 5 quantifies deviation from the standard self similar model (see Equat ion 



(9|). We give in Table |EiJ (|E2 I 



0. Similarly, Table 



E3 



E8 we show 



and Cp in the case of all clusters (resolved clusters), with S 
(E4l shows our measurements in the case of all clusters (resolved clusters), with 5 = 0.12. In Tables E5 
the result of our analysis on mass bins 1 — 4 introduced in Table [Ta[ respectively (with 3 = 0). Similarly, Tables [E9| - )E11| 
provide the best fit universal pressure values and the corresponding covariance matrix for mass bins 1 — 3 introduced 



in Table lb respectively (with 6 = 0). 
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Figure Dl. Spectrum of tiie covariance matrix and the contribution of different eigenmodes to Xq- 
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Figure D2. Components T„ j; of the eigenvectors Tn of the covariance matrix. Here n labels different eigenvectors, chosen so that 
the eigenvalues satisfy Ai < A2 < ■ • • < Ajv^. The label fc runs over all radial bins. 

Table El. Pressure measurement of all MCXC clusters with the standard self-similar model (<5 = in Equation l[9|l). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 


Pressure 


2904.845 


503.878 111.528 -8.831 8.515 


54.610 


-21.088 8.990 






(a) Best fit pressure values P (xlO" 


■^) 






Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 



Bin 1 


359.047 


-82.359 


16.095 


-2.989 


0.872 


-0.540 


0.220 


-0.016 


Bin 2 


-82.359 


36.703 


-12.792 


2.632 


-0.470 


0.052 


-0.023 


0.018 


Bin 3 


16.095 


-12.792 


10.142 


-4.167 


0.720 


-0.010 


-0.067 


0.031 


Bin 4 


-2.989 


2.632 


-4.167 


4.320 


-1.928 


0.298 


0.010 


-0.017 


Bin 5 


0.872 


-0.470 


0.720 


-1.928 


2.405 


-1.162 


0.200 


-0.008 


Bin 6 


-0.540 


0.052 


-0.010 


0.298 


-1.162 


1.577 


-0.786 


0.132 


Bin 7 


0.220 


-0.023 


-0.067 


0.010 


0.200 


-0.786 


1.094 


-0.476 


Bin 8 


-0.016 


0.018 


0.031 


-0.017 


-0.008 


0.132 


-0.476 


0.432 



(b) Covariance matrix Cp (xlO ^) 
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Table E2. Pressure measurement of resolved MCXC clusters with the standard self-similar model ((5 = in Equation Jol). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 


Pressure 


3156.051 


652.465 163.366 2.700 48.387 


78.075 


-14.598 14.782 






(a) Best fit pressure values IP (xlO" 


-') 






Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 



Bin 1 


385.004 


-76.888 


18.073 


-2.839 


1.073 


-0.664 


-0.018 


-0.199 


Bin 2 


-76.888 


42.837 


-11.045 


3.328 


-0.439 


-0.052 


-0.177 


-0.240 


Bin 3 


18.073 


-11.045 


11.960 


-3.454 


0.846 


-0.087 


-0.154 


-0.136 


Bin 4 


-2.839 


3.328 


-3.454 


5.082 


-1.641 


0.291 


-0.054 


-0.141 


Bin 5 


1.073 


-0.439 


0.846 


-1.641 


2.785 


-1.022 


0.214 


-0.084 


Bin 6 


-0.664 


-0.052 


-0.087 


0.291 


-1.022 


1.812 


-0.677 


0.102 


Bin 7 


-0.018 


-0.177 


-0.154 


-0.054 


0.214 


-0.677 


1.281 


-0.466 


Bin 8 


-0.199 


-0.240 


-0.136 


-0.141 


-0.084 


0.102 


-0.466 


0.745 



(b) Covariance matrix Cp (xlO ^) 
Table E3. Pressure measurement of all MCXC clusters with the modified self-similar model (5 = 0.12 in Equation Kf). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 


Pressure 


2819.487 


515.355 97.656 -15.799 4.248 


50.901 


-22.110 8.053 






(a) Best fit pressure values IP (xlO' 


"') 






Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 Bin 8 



Bin 1 


352.226 


-82.197 


16.578 


-3.117 


0.820 


-0.496 


0.211 


-0.018 


Bin 2 


-82.197 


36.333 


-12.687 


2.618 


-0.467 


0.035 


-0.018 


0.023 


Bin 3 


16.578 


-12.687 


9.842 


-4.026 


0.683 


0.001 


-0.075 


0.036 


Bin 4 


-3.117 


2.618 


-4.026 


4.114 


-1.830 


0.269 


0.019 


-0.019 


Bin 5 


0.821 


-0.467 


0.683 


-1.830 


2.273 


-1.093 


0.179 


-0.004 


Bin 6 


-0.496 


0.035 


0.001 


0.268 


-1.093 


1.487 


-0.740 


0.121 


Bin 7 


0.211 


-0.018 


-0.075 


0.019 


0.179 


-0.740 


1.034 


-0.451 


Bin 8 


-0.018 


0.023 


0.036 


-0.019 


-0.004 


0.121 


-0.451 


0.405 



(b) Covarianee matrix Cp (xlO '^) 
Table E4. Pressure measurement of resolved MCXC clusters with the modified self-similar model {S = 0.12 in Equation Kt). 





Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 


Pressure 


3095.361 


672.793 


159.082 


0.880 


46.378 


75.526 


-16.262 8.843 






(a) Best fit pressure values IP (xlO" 


-') 








Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 


Bin 1 


378.015 


-77.464 


18.369 


-2.951 


1.042 


-0.619 


-0.015 


-0.185 


Bin 2 


-77.463 


42.440 


-11.029 


3.331 


-0.412 


-0.058 


-0.166 


-0.237 


Bin 3 


18.369 


-11.029 


11.612 


-3.340 


0.811 


-0.060 


-0.153 


-0.126 


Bin 4 


-2.951 


3.331 


-3.340 


4.854 


-1.554 


0.263 


-0.036 


-0.140 


Bin 5 


1.042 


-0.412 


0.811 


-1.554 


2.639 


-0.959 


0.193 


-0.072 


Bin 6 


-0.619 


-0.058 


-0.060 


0.263 


-0.959 


1.713 


-0.638 


0.094 


Bin 7 


-0.015 


-0.166 


-0.153 


-0.036 


0.193 


-0.638 


1.214 


-0.442 


Bin 8 


-0.185 


-0.237 


-0.126 


-0.140 


-0.072 


0.094 


-0.442 


0.711 



(b) Covariance matrix Cp (xlO '^) 
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Table E5. Pressure measurement of Mass-bin 1 of all MCXC clusters, as defined in Table [la|(<5 = in Equation Kf). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 


Bin 8 


Pressure 


2864.897 


143.372 165.443 82.563 57.910 


-22.757 


53.503 


37.140 






(a) Best fit pressure values P (xlO" 


-^) 








Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 


Bin 6 


Bin 7 


Bin 8 



Bin 1 


1351.387 


-266.076 


30.865 


0.085 


1.519 


-0.727 


-0.089 


0.199 


Bin 2 


-266.075 


117.062 


-37.249 


5.581 


0.115 


0.093 


-0.088 


-0.025 


Bin 3 


30.866 


-37.249 


34.247 


-13.720 


2.080 


0.136 


0.038 


-0.086 


Bin 4 


0.084 


5.582 


-13.720 


16.293 


-7.307 


1.239 


0.070 


-0.055 


Bin 5 


1.520 


0.115 


2.080 


-7.307 


9.706 


-4.636 


0.856 


-0.025 


Bin 6 


-0.728 


0.093 


0.136 


1.239 


-4.636 


6.452 


-3.190 


0.559 


Bin 7 


-0.089 


-0.088 


0.038 


0.070 


0.856 


-3.190 


4.424 


-1.869 


Bin 8 


0.199 


-0.025 


-0.086 


-0.055 


-0.025 


0.559 


-1.869 


1.832 



(b) Covariance matrix Cp (xlO ^) 



Table E6. Pressure measurement of Mass-bin 2 of all MCXC clusters, as defined in Table [Ta|(<5 = in Equation Ky). 





Bin 1 


Bin 2 Bin 3 Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 


Pressure 


2546.299 


226.786 253.102 18.308 


-79.667 


204.195 


-123.737 


-7.187 






(a) Best fit pressure values P (xlO" 


-') 








Bin 1 


Bin 2 Bin 3 Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 



Bin 1 


2504.452 


-719.122 


230.700 


-84.783 


37.468 


-19.648 


10.328 


-3.445 


Bin 2 


-719.123 


311.542 


-137.925 


54.001 


-22.478 


10.146 


-4.269 


1.271 


Bin 3 


230.702 


-137.925 


99.347 


-53.731 


22.694 


-9.263 


3.633 


-0.946 


Bin 4 


-84.785 


54.002 


-53.731 


46.988 


-27.135 


11.806 


-4.740 


1.355 


Bin 5 


37.469 


-22.479 


22.695 


-27.135 


26.089 


-15.997 


6.810 


-1.871 


Bin 6 


-19.648 


10.146 


-9.264 


11.806 


-15.997 


16.275 


-9.554 


2.819 


Bin 7 


10.328 


-4.269 


3.633 


-4.740 


6.809 


-9.554 


9.325 


-3.905 


Bin 8 


-3.445 


1.271 


-0.946 


1.355 


-1.871 


2.819 


-3.905 


2.673 



(b) Covariance matrix Cp (xlO ^) 



Table E7. Pressure measurement of Mass- bin 3 of all MCXC clusters, as defined in Table [Ta| (i5 = in Equation Ky). 





Bin 1 


Bin 2 Bin 3 Bin 4 


Bin 5 


Bin 6 


Bin 7 Bin 8 


Pressure 


2614.687 


599.174 -119.300 -23.585 


9.475 


161.234 


-148.958 139.492 






(a) Best fit pressure values P (xlO" 


-') 






Bin 1 


Bin 2 Bin 3 Bin 4 


Bin 5 


Bin 6 


Bin 7 Bin 8 



Bin 1 


816.335 


-142.117 


21.841 


-11.060 


7.061 


-4.547 


1.150 


0.328 


Bin 2 


-142.117 


90.898 


-31.204 


6.198 


-3.173 


1.336 


-0.726 


0.394 


Bin 3 


21.841 


-31.204 


33.339 


-13.358 


2.133 


-0.887 


0.478 


-0.212 


Bin 4 


-11.060 


6.198 


-13.358 


16.314 


-6.910 


1.096 


-0.559 


0.262 


Bin 5 


7.061 


-3.173 


2.133 


-6.911 


9.575 


-4.265 


0.700 


-0.205 


Bin 6 


-4.547 


1.336 


-0.887 


1.096 


-4.265 


6.279 


-2.804 


0.305 


Bin 7 


1.150 


-0.726 


0.478 


-0.559 


0.700 


-2.804 


4.368 


-1.827 


Bin 8 


0.328 


0.394 


-0.212 


0.262 


-0.205 


0.305 


-1.827 


1.997 



(b) Covariance matrix Cp (xlO ^) 
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Table E8. Pressure measurement of Mass-bin 4 of all MCXC clusters, as defined in Table [la|(i5 = in Equation dOl). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 


Pressure 


4522.423 


811.544 140.014 -63.051 -16.867 10.731 


60.649 


-94.528 


(a) Best fit pressure values IP (xlO^^) 




Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 



Bin 1 


4322.774 


-1224.992 


307.886 


-40.992 


-10.662 


8.444 


-1.387 


-0.802 


Bin 2 


-1224.992 


434.752 


-126.640 


19.097 


2.756 


-3.787 


0.775 


0.135 


Bin 3 


307.888 


-126.641 


59.055 


-15.870 


0.103 


1.643 


-0.966 


0.207 


Bin 4 


-40.993 


19.097 


-15.870 


13.835 


-4.773 


-0.305 


0.650 


-0.317 


Bin 5 


-10.663 


2.756 


0.103 


-4.773 


6.955 


-2.669 


-0.122 


0.219 


Bin 6 


8.445 


-3.788 


1.643 


-0.305 


-2.669 


4.461 


-1.854 


0.077 


Bin 7 


-1.387 


0.775 


-0.966 


0.650 


-0.122 


-1.854 


3.298 


-1.417 


Bin 8 


-0.803 


0.135 


0.207 


-0.316 


0.219 


0.077 


-1.417 


1.642 



(b) Covariance matrix Cp (xlO ^) 



Table E9. Pressure measurement of Mass-bin 1 of resolved MCXC clusters, as defined in Tabic [Tb|(<5 = in Equation Kf). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 


Pressure 


2970.143 


-124.936 58.889 102.331 97.592 54.978 


75.993 


94.229 


(a) Best fit pressure values IP (xlO~^) 




Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 



Bin 1 1375.649 -236.216 47.033 



Bin 2 
Bin 3 
Bin 4 
Bin 



-236.214 129.290 



-31.665 



-2.420 
8.192 



2.119 



-1.271 



-1.240 -0.716 



-0.684 -0.393 -0.512 -0.764 



5 
in 6 

in 7 
in 8 



47.036 
-2.420 
2.120 
-1.271 
-1.240 
-0.716 



-31.667 38.119 



-11.362 2.634 



-0.653 -0.387 -0.557 



8.192 
-0.685 
-0.392 
-0.512 
-0.764 



-11.361 
2.634 
-0.653 
-0.387 
-0.557 



17.521 
-6.171 
1.401 
-0.481 
-0.385 



-6.171 
10.149 
-4.066 
1.023 
-0.461 



1.401 
-4.066 
6.674 
-2.656 
0.468 



-0.481 
1.023 
-2.656 
4.630 
-1.700 



-0.385 
-0.461 
0.468 
-1.700 
2.578 



(b) Covariance matrix Cp (xlO '^) 



Table ElO. Pressure measurement of Mass-bin 2 of resolved MCXC clusters, as defined in Table [Tb| {S = in Equation l|9Jl). 





Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 


Pressure 


2482.111 


747.796 


-29.912 


-32.605 


2.269 


180.980 


-115. 17C 


74.383 






(a) Best fit pressure values IP (xlO" 


-') 








Bin 1 


Bin 2 


Bin 3 


Bin 4 


Bin 5 


Bin 6 


Bin 7 


Bin 8 


Bin 1 


682.399 


-103.800 


28.794 


-8.645 


4.760 


-4.283 


0.648 


-0.521 


Bin 2 


-103.799 


88.006 


-19.527 


7.560 


-3.040 


0.400 


-1.156 


-0.276 


Bin 3 


28.794 


-19.527 


31.489 


-8.677 


2.714 


-1.348 


-0.072 


-0.722 


Bin 4 


-8.645 


7.560 


-8.676 


15.190 


-4.710 


1.212 


-0.803 


-0.205 


Bin 5 


4.760 


-3.040 


2.714 


-4.710 


8.745 


-3.022 


0.850 


-0.533 


Bin 6 


-4.283 


0.400 


-1.348 


1.212 


-3.022 


5.652 


-1.930 


0.398 


Bin 7 


0.648 


-1.156 


-0.072 


-0.803 


0.850 


-1.930 


3.906 


-1.247 


Bins 


-0.521 


-0.276 


-0.722 


-0.205 


-0.533 


0.398 


-1.247 


2.349 



(b) Covariance matrix Cp (xlO '^) 
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Optimal SZ measurement 25 

Table Ell. Pressure measurement of Mass-bin 3 of resolved MCXC clusters, as defined in Tabic [lb] (5 = in Equation k\). 





Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 


Pressure 


3737.228 


1220.556 334.228 7.300 34.974 23.814 


7.113 


-101.074 


(a) Best fit pressure values P (xlO^^) 




Bin 1 


Bin 2 Bin 3 Bin 4 Bin 5 Bin 6 


Bin 7 


Bin 8 



Bin 1 3842.297 -1117.928 295.139 -49.881 -6.830 7.859 -1.586 -0.158 

Bin 2 -1117.928 417.687 -116.036 26.214 2.461 -3.534 0.417 -1.168 

Bin 3 295.139 -116.036 58.291 -14.939 1.537 1.550 -1.007 -0.033 

Bin 4 -49.882 26.214 -14.939 14.707 -4.168 0.113 0.563 -0.664 

Bin 5 -6.829 2.461 1.537 -4.168 6.989 -2.360 0.175 0.093 

Bin 6 7.858 -3.534 1.550 0.113 -2.360 4.449 -1.665 0.203 

Bin 7 -1.586 0.417 -1.007 0.563 0.175 -1.665 3.262 -1.232 

Bin 8 -0.158 -1.168 -0.033 -0.664 0.093 0.203 -1.232 2.008 

(b) Covariance matrix Cp (xlO~^) 
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